STOKES MATRIX APPROACH TO 
CLUTTER MODELLING 


A Thesis Submitted 

In Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


if ft 8 

' * < ' ' ~ \ A * t 

by * - " 

P N SRIKANTH 


to the fa 

DEPARTMENT OF ELECTRICAL ENGINEERING 
1M>I \N |\V| I'll 'll 01 IL(JIV)I(K;Y. K\NI»1JK 

SEPTEMBER. 1984 


/ 


A V 

\ 


S. 


V 



CERTIFICATE 


Certified that the work reported m the 
thesis titled "STOKES MATRIX APPROACH TO 
CLUTTER MODELLING" has been carried out under 
our supervision by Mr P.N.Srikanth and has 
not been subm tted elsewhere for a degree. 



Dr P.R.K.Rao 


Professor 

Dept of Elect Engg 
IIT Kanpur 


Dr K.R.Sarma 
Professor 

Dept of Elect Engg* 
IIT Kanpur 


> 


POT GRAPH \ I F, OFFICE 
I Ins (bests h s b > n ‘pptoyed 
for the m ar i <1 lb 1 Xf s'tec of 
\I i f I u 1 ill !htOi< .v ( t 1 COCl ) 


in «l ecu 


la 1 la Vu.lt t '»«, 


rrguUtions jf the Indian 
Institute of lechuologv 


D.tfd 


m $ 




tanprtr 

Li 



i\ m* 


i ) 4 
f 4 mi' 

*■ ** * 


^ * 

% 



EE-t m-^~^~ S7D 


i 


JA h n ) 1 

H I*t/ii A ,iMl tiVH ill i 

wjtr nr • 



ACKNOWLEDGEMENTS 


I am greatly indebted to my guides Dr. K.R. Sarma and 
Dr. P.R.K. Rao for their invaluable guidance and constant 
encouragement extended to me during my thesis work. 

I also thank Dr. N.C. Mathur and Dr. V.P. Smha for th< 
help m solving some of the problems I faced m this work. 

I thank all of my friends who made my stay at I,I*T* 
Kanpur enjoyable. 


- P.N. SRIKANTH 



ABSTRACT 


This thesis is an attempt to model the electromagnetic 
returns from clutter-producing targets that degrade the 
Radar detection performance. 

The study is carried out In the frame-work of scattering 
theories employing target and Stokes reflection matrices. 

The use of Stokes reflection matrix allows for the decompo- 
sition of a target into a deterministic target and a noise 
target. The deterministic part of the target model admits 
a characterization of a target m terms of five phenomeno- 
logical parameters orientation, helicity, skip and character- 
stic angles and 'magnitude' of the target. 

Expressions for target scattering and Stokes reflection 
matrices for static surfaces modelled as random rough surfaces 
are derived. These are extended to obtain expressions for 
doppler spectrum of the scattered electromagnetic returns when 
temporal variations of the surfaces are present. 

To illustrate these ideas, an asphalt road and a sea- 
surface are taken as examples. The variation of five pheno- 
menological parameters associated with above two targets with 
frequency, aspect angle and surface roughness are also 


studied 
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CHAPTER 1 
JiiT RODUCTIOn 

1.1 motivation and background 

The received signal m radar due to scattering 
of electromagnetic waves from ob]ects othei than the 
desired one is termed as clutter. For better perfor- 
mance of the raaai receiver, the effecls due to clutter 
must be minimized. In order to minimize these effects, 
one must know the characters tics of the clutter. The 
experimental data about the clutter characterstics 
allow one to design radar receivers. But these data 
are non-stationary and varies with geographical and 
radar parameters. Hence it is desirable to have a 
theoritical model for clutter. This gives insight into 
the clutter phenomena, m addition to the design of 
adaptive radar recei/ers. 

A number of clutter models are available m 
literature. One such model considers the target ae a 
collection of ellipsoids with varying cross-sections[5] . 
”|he scattered power return due to a single ellipsoid 
was calculated first and the total scattered power is 
computed as the some of all the contributions from 
individual ellipsoids. But such an approach is fundamenta 



incorrect, as the presence of other ellipsoids will 
effect the scatterang characterstics of individual 
ellipsoids. Wong et.al [4] modelled the clutter as a 
random collection of rotating dipole scatterers* The 
rotation of dipoles represents the random motion of 
the objects such as leaves, vegetation, branches etc. 
due to wind forces. But the expressions obtained using 
this model for doppler spectrum is not related to any 
of the physical characterstics of the surface like its 
shape, roughness or permittivity. It is normal practice 
m detection and estimation theory to assume some kind of 
clutter distribution like Gaussian, log-normal or weibull 
for the design of Radar receivers [6]. In these, some 
parameters are introduced whose values are obtained by 
fittm f with experimental data. But these parameters are 
also not related to any of the properties of the surface. 
Hence it is desirable to obtain a model which contains 
parameters which are related to the physical properties 
of the surface. 

The characters tics of clutter can also be studied 
from its target scattering matrix (TSk) and Stokes 
reflection matrix associated with it. These two matrices, 
express the scattering and polarization properties of 
the target. Thus the information content m the pola- 
rization which was neglected m the past is utilized 
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to study the clutter characterstics. The study of TSM 
and Stokes reflection matrix reveals some interesting 
properties of the target. The eigenvalues and eigen- 
vectors associated with the TIL is related to some of 
the properties of TSM through five phenomenological 
paramotcis proposed by Huyncn m. The basic idea is 
to decomposed a target return into one due to a deter- 
ministic part of the targ >t and a noisy part. The TSM 
corresponding to the deterministic target is obtained 
and the associated five parameters represent clutter 
characterstics like its symmetry, non-symmetry, size 
and orientation. 

TSiv and Stokes reflection matrices can be obtained 
exoerimen tally. Huynen has suggested an experimental 
procedure to measure the TSjvil" 3 ] using duaJ polarized 
antennas, for details of experimental procedure, we 
reftr to the aoove cited papor-by Huynen. Theoretically, 
these matrices can be obtained theoretically utilizing 
scattering tneories. Scattering theories are well 
developed m the literature. In fact, several excellent 
books are devoted to this topic. This includes the book 
by Beckmann and Spizzichmo [7] and Bass and Fuks [9]. 
These scattering theories are utilized to model clutter 
especially sea clutter[l7]. Apart from clutter modelling 
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scattering theories also find application m the area 
of remote sensing [22] . 

In this thesis, we compute the scattering 
matrix and Stokes reflection matrix associated with 
a surface modelled as random rough surface. The 
Stokes reflection matrix thus obtained is decomposed 
into one due to the deterministic part of the target 
and other due to the noisy part of the target. The 
TSm associated with deterministic part is obtained 
and the five phenomenological parameters are cal- 
culated from TSM, 
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1.2 ORGANIZATION OF THE THESIS : 

In chapter 2, we review some of the concepts of 
polarizanon and its repres -ntation . The properties 
of TSJv and otokcs refine tion matrices are discussed m 
detail. 


In chapter 3, we obtain expressions for mean and 
variance of TSJu for a static surface modelled as a 
random ' tou§n surface. These expressions 

are then extended t c include the temporal variations 
of the surface. 

In chapter 4, we numerically compute the Stokes 
reflection matrices for an asphalt road and sea-surface. 
The TSJvi corresponding to deterministic part of the 
target is obtained and the associated target para- 
meters are computed. The variation of these para- 
meters with aspect angle, a frequency and surface rough- 
ness are also studied. 

The results obtained m chapters 2,3, and 4 
are discussed and concluded upon m chapter 5. 
Suggestions are given for the further improvements 
of the present model* 





CHAPTER 2 
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POLARIZATION TARGET SCATTERING AND 
STOKES REFLECTION MATRICES 


In this Chapter we discuss the concept of polarization 
and its representation. This is followed by the dia&ssion 
of some properties of target scattering matrices (TSM) and 
their relationship with five phenomenological parameters 
and stokes ref lection matrices . 

2.1 Polarization and its representation 

Let E and H denote the electric and magnetic fields 
of a electromagnetic (e.m.) wave at a point m a medium with 
associated poyntmg vector S m the Z-directon of the 
rectangular co-ordinate system XYZ. Then the E and H fields 
lie m the transversal plane x-y. Let x and y components 
of electric field E be 

E x = E xo Cos (wt +S 1 > 7 E y = E yo Cos(wt+€> 2 > ...(2.1.1) 

where =w^,-kz and G — u -kz with 

k describing the wavenumber associated with the medium 
Denoting the phase-difference ( ) bet ' E x and E^. 

components by & , the following can be shown to hold good. 



Ex' 


Ex 2 o 


Cos'S 

xo yo 


+ 


E 2 

E 2 

yo 


= sun 


2 S 


or 


aE^ - bE x E y + CEy 


= 1 


where a = 


E xo sin 




; b = 2 c OS § , o= 1 

E xo E yo sin ^ El^sm 2 S 


yo 


The ellipse specified by eqn (2.1.2) is called 
polarization ellipse and is described by following three 
parameters . 

i ) ellipticity angle X = tar? ( - 9— - ) f - 45 ?, "X S. 45°J 

ii) magnitude of the ellipse m = ^ OA 2 +OB 2ft M ___ (2 .1. 3 ) 

iii) orientation angle of the ellipse = vp 

v. (see Fig. 2.1.1.) 

These three parameters Vn , \y and X completely describe 
the polarization ellipse and are related to E^, e 2 , S 
as under ; 


m 


=c e i + 


tan 2*i>- 


sm 2T. = 


2E 1 E 2 


E 1~ E 2 

2E 1 E 2 


CO 


sS 


...(2.1.4) 


_ 2 „2 
E l +E 2 


sin£ 
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2 

It may be noted that m 
power m the e.m. wave. 


is a measure of the total 


For the purposes of this study, instead of the repre- 
sentation for the electric field given m (2.1.1), an 
alternative representation m terms of m. and is more 
convenient.^ e develop this alternative representation for 
the case of zero orientation angle first and then for arbitary 
angle 


For § =%/2 
’e 


E 


L E y 


, we have V 
j~Ei cos( w "t+9^) 

E 2 cos( w "t + ®2^ 


= o 


L 


Re 


[E 1 e ^ w 't+© 1 ) 
e (j wt+© 2 ) 


= Re 


l 1 


E 2 e 


1 $ 




Ra 


lEi 


jE 2 


e ^3 (ot+9^) 


Substituting E 1 = mcost, ; E£ =-msinT , we get in 
complex notation. 


E = m 


Cosx 


10 

e JW i 


-j smX 

where following usual practice we have supressed exp 
(jwt) term. 


The representation for a arbitary orientation angle 
corresponds to a rotation of the x-y co-ordinates thro' an 
angle ty/ . ^nd E may be seen to be given by 
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E (a,V<,-C ) = 


(cosy - sin^l 


'cos X * 

m I 

tsmy cosy 


jsmx 


>J 9 1 


where 


is the absolute phase. 


..(2.1.5) 


The above representation is often called the geometrical 
representation . 


Some of the special cases of interest are - 

If E^O, the wave is linearly polarized in x - direction. 

If E^o, the wave is linearly polarized in y - direction 

O <r 

If 0 - ^ wave 1S linearly polarized with 

45° orientation. 

«. and E^=E 2 

If &=±90° i /, the wave is circularly polarized. 

4 = +90° represents left circularly polarized wave and 
§ = —90° represents right circularly polarized wave. It may 
be noted that our convention about left and right circular 
polarizations are in accordance with IEEE standards is at 
variance with the practice m conventional optics. 

The geometrical parameters and t can be represented 

on a sphere called Poincare sphere [ 1 i . In this sphere 
each polarization state is represented by a unique point. 

Tne size of the sphere depends on the magnitude of the 
ellipse and the latitude and longitude of a point represent 
and 2X. associated with E at that point. The following 



situations are of special interest 


11 


(1 ) The equator of the ©oincare sphere represents 
all linearly polarized waves (2T = 0°) 

(2) The poles represent circularly polarized waves 
(2y= jfc 90°) . 

(3) The upper hemisphere represents all left sensed 
polarizations and lower hemisphere represents all 
right sensed polarizations. 


( 


We may also associate stokes vector with polar izationL |S Q 
) . The stokes vector tj is defined as 


2 = (I, Q, U, V)' 


where 


I = e; 


+ 



Z - S 


U = 2 E 1 E 2 cos^ 


( 2 . 1 . 6 ) 


V = 2 _ E.E.sin^ 

z 

mtrisi ic 

where Z is the 'lmpedence of the medium. 

In terms of geometrical parameters, these components of g 
are given by 


I 


= : Q = m^cos2<#/ cos2T 

U = m^sw2^ cos2T. ; V = m^sin^2t 


i2 

(2.1.7) 


From the above relations, it can be seen that, 

2 2 2 2 
I = Q +lT=fV 


( 2 . 1 . 8 ) 


Thus only 3 of the 4 parameters are independent. 

An electromagnetic wave satisfying the above condition is 
said to be completely polarized. An e.m. wave which 
does not satisfy above condition is said to be partially 
polarized or unpolarized. 


-In above discussions, E^E^ are taken as constants. 
Let us now take them as random samples of a stationary, 
erogodic, random process. In this situation, the equatms 
(2.1.1. - 2.1.3) are modmed as follows 

Ex = E (t) cos (wt-ta (t)); 

XO Jl 

E = E (t) cos (wt+0i (t) ) + & (t)) 
yo yo 1 


and the equation of ellipse becomes, and the equation 


of the ellipse becomes, 

A(t)Ex - b(t)E x E y + c(t)Ey = 1 


..( 2 . 1 . 8 ) 


where a(t) = 


E?(t) 


sin 2 ^ (t) 
c(t) = 


b(t) ~ - l- cos ? 

E^ (t )E 2 (t)sin^(t) 


(t)sin z 4 (t) 
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Thus the major and minor axis in this case are time varying. 
It corresponds to a wave whose orientation, eccentricity and 
size of the polarization ellipse continuously changes with 
time. The Stokes parameters for this case are given by 
the following time-averages ' > , 


I = 


Q = 


U = 

V = 

where ^ 


<E^(t)> <E^(t)> 

— + — 

Z Z 

< E^ (t)> < E^CtV 

Z Z 

^ <.E^(t)E 2 (t) cos 3 (t)> 
rjr v E^tJE^lt) sm S (t)> 
> indicates time averaging. 


... (2.1.9) 


It may be shown that - 
2 2 2 2 

I Q +U +V and e.m. waves satisfying this 

inequality strictly are called partially polarized provided 
all the quantities Q,U,V are not simultaneously equal to zero. 

For unpolarized waves, E^(t) and B^it) are uncorrelated and 
hence U=V*o. Since Q =m cos 2^cos 2C and average value 
of </cos 24 /cos 2 O = 0, Q=0, Hence £N=U=V=0 for a unpolarized 


wave 
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Thus the 0,11, V components represent the polarized 
partof the wave. If anyone of them is non-zero, it indicates 
that the wave contains a polarized component. 


The degree of polarization is defined as the ratio 
of polarized power to total power. Thus 


degree of polarization = P°l;rized power 

tbtal power 


i .e , 


2 2 2 u 

d = (Q tlT+V^ 

I 


( 2 . 1 . 10 ) 


The auto correlation and cross correlation 
of E-^Ct) and @ 2 (t) are often used to define the complex 
correlation matrix (a) so called coherence matrix J 03 
i as follows: 


J = 


XX 


xy 


< E 1 (t)E* (t) >< E 1 (t)E*(t) > 


— ( 2 . 1 . 11 ) 


yiz 


yy 


<E 1 (t)E 2 (t)> <E 2 (t)E*(t) > 


The complex degree of coherence Ury is defined as follows 


^ xy = < 


J *y 


( xx ) 2 ( yy ) 2 


= (^xy) 


© 


xy 


.. ( 2.1 . 12 ) 


where ^xy & 1 

when M'xy = 1 , the wave is said to be completely polarized, j 
when f$xy 4k 1/ the wave is said to be partially polarized, j 
when ^xy = 0, the wave is said to be unpolarized. 
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Thus the concept of temporal coherence and polari- 
zation are intimately connected. 

Any partially polarized wave with Stokes vector 

jp 

['I/Q,U,V] may be considered to consist - of an unpolanzed 

2 2 2 2 2 

wave with an intensity I =1- (Q+U+V ) and a completely 

o 2 2 2 

polarized wave with an intensity I_=Q+U+V. Thus 
2 2 2 

I = I u + I is the sum of intensities of the polarized and 

2 2 2 2 

unpolarized waves, and I ^ Q + U + V with equality sign 
holding good for completely polarized wave. This kind of 
decomposition of a partially polarized wave into a polarized 
wave and an unpolarized wave is utilized in the target 
decomposition theorems derived by Huyhen [ i ] . 

We conclude tus discussion by listing the stokes 
parameters of some commonly encountered polarizations 

1 linearly polarized wave with E 2 =0 i»e.# 

I = P, Q=P; U=V=0 

2 linearly po'arized wave with E^=0; (le) 

I = P, Q=P, U=V=0 

3 left circularly polarized wave with E 1 =E 2 **90° 

I = P; Q=0, U=0, V=P, 

Where m a h ove eqns P denotes power density. 


2.2 Target Scattering Matrix 

Let us consider an e.n. wave E 1 incident 
on a point target. The scattered wave E s from 
that point is given by 

E S = SE 1 (2.2.1) 

where S is called scattering matrix and is a 2x2 
complex matrix. 

( 2 . 2 . 2 ) 

where s^ 1 = E® I 



In our work, we shall consider only targets with 2x2 
complex symmetric matnc's (i.e.) s^ = s 2 i 

In our work, we consider the orthognal, Honzental 
and vertical Polarizations pair (H-V) to define the 
scattering matrix S. This matrix S can easily be converted 
to any other matrix associated with some other orthogonal 
polarization pair base by simple transformations^! fl . 


Let S = 


11 

s 12 

21 

S 22 
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The polarization scattering properties of a 
target are contained m its TSM. The properties of this 
matrix, can be investigated conveniently using its 
eigenvalues and eigen vec Lors , Let us formulate eigen- 
value problem as follows 0 

S x = t / ... (2.2.5) 

Thfs equation has two solutions. Denoting by t^,t 2 

eigenvalues and by xl, x_2, the eigenvectors of (2.2.5) 
we have, 

S _x_l = and S x2 = t 2 jc2*_ ... (2.2.6) 

The emgen vectors j<l and _x2 are orthogonal and are 
normalized. 

The method of computing these eigenvalues and eigen- 
vectors are sligntly different from computing the con- 
ventional eigenvalues and eigenvectors. In general, 
direct solution of eqn. (2,2.5) is difficult. Hence we 
have to compute these values m a slightly different way, 
which we vill discuss m 

The scattered field ~ is given 

E s = SE 1 where E 1 is the incident field. 

The scattered power Ps is .jiven by - 


by - 
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Ps = Ef + E S = (SE 1 ) "^SE 1 ) ...(2.2.7) 

expressing E 1 mterms of xl and x2 , 

E 1 = E^ Xl + E 2 *2 ...(2.2.8) 

substitution in eqn. (2.2.8) gives, 

PS = ( SE^ XI + SE^ x2) + (SE^ Xl + SE| x2) 

knowing that SXl = txl* and S X2 = tx2* and 
using the property that xl and x2 are orthonormal vectors, 
pwe get, 

P s = Kl 2 H 2 + E 2 P 2 2 

P A Wl ~ | B 2 j 2 (jt^ 2 -^t 2 | 2 ) ...(2.2.9) 

where p^ - | + | E 2| = 91 A = total tiansunt^ power 

Without loss of generality, we assume I 'z l^l 
and let be fixed. 

It can be easily seen that eqn. (2.2.9) is maximum when 

1 E |i - 0 

Tli is Pma^v == | t^ \ 2 . The value of | t^l is 

assigned as a parameter . m which is called target magnitude. 

It may be viewed as a overall measure of the target size. 

From above discussions, it can be seen that the 
principal eigen vector xl corresponds to polarization 
for which the return power is maximum. xl can be represented 
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in terms of geometrical parameters associated with it 
as follows - 


The V^and associated with this max. polarization is 

assigned to the target and are called target orientation 
and target helicity respectively , The eigen values are 
m general, complex and can be expressed in the following 
form -- 

t ^ = m e J ®l and t 2 = am .. (2.2.12) 

where the conditions j J *= m and jt^ j ^ j t 2 j if 
a 41 are satisfied. 

Having found eigenvalues and eigenvectors, let us 
see about nulls associated with T.S.M. There are two hinds 
of nulls - Co-POL nulls and X-POL nulls. They are defined 
as follows 

The Co— POL nulls is polarization of radar n of 
identical Tx and Rx antennas which produces zero voltage 
reception at: receiver terminals. 

null 

The X-POL^/is that polarization n, of Orthogonal 
TX and RX antennas, which produces zero voltage tece»*ption 
at the receiver terminals. 
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Mathematically, n' sn =0 for CO-POL nulls 

nj_sn =0 for X-POL nulls *♦. (2.2.1 

where n, n^are orthogonal polarizations. 

CO-POL and X-POL nulls are practically measurable 
quantities. Hence it is dcsireable to express these nulls 
m terms of eigen values and eigenvectors of T.S.M. 
Expressing polarization n interns of X 1 & x 2. 


Let n = njjcl + n ^ x2 


.. . (2 .2.13) 


Substituting m equation for Co- POL nulls, 

(n 1 xi + n 2 x2 ) (Sn 1>c i + Sn 2 x2) = 0 

using Sx = tjj* and orthonormal properties of X 1 and 
X 2 , one finally gets. 


£l = - 2 (a)^ e ^ (e 2-e,.)/2 £i + 22 

n 2 = + jfa)* 5 e j( ®2- e l )/2 21 + x 2 


... (2.2 .14) 


51 

5 2 


expressing m symmetrical form - 

= (a)% e ~ j45 ° e 3 (9 2^ 6 1 )//4 x ± + e j45 ° 

p _ 45 0 e J (®2~®1 '^ //4 x + t 45° 

= -(o.) 2 e J4b 6 -1 e 3q ^ 


e -j ( © 2 -© 1 )/4^ 
e -j -Q x )/4 X 


At this stage, for convenience, let us say oc 


tan 2 Y 


and 



Thus and n 2 


21 


n 2 


becomes. 


-#(45°+V). 


= Bin y e* ' 'x^ + cos Ye"’" x 


j (45°+», 


-2 


= . S m y e -j (45 u +‘^)x 1 + cosy e J (45 +>>) x 


(2.2.15) 


where y is called character Stic angle. As can be seen, the 
angle subtended by ?nl , n2 at the centre of the pomcare’ 
sphere is 4 y. Hence parameter y dahas acquired physical 
significance. The parameter is called "target skip 
angle" and we will discuss the role of this parameter later. 

The X-POLS nulls can be easily seen to be eigen- 
vectors Xl and X2 themselves. 


It is interesting to note that the polarizations nl, n2, 
xl, x2 when plotted in the foincare 1 sphere lie on one 
great circle. The line joining Xl and x2 bisects the 
angle subtended by nl and n2, when these points are connected 
to the c entre of the circle, they form the shape of a fork 
which is called Huynen's fork. From eqns (2.2.15), one 
can see that the parameter determines the rotation the 
* prongs’ of the fork about line joining xl, x2 by an 
amount of 2%) 

Given the CO-POL and X-POL nulls, one can determine 
TSM except for target magnitude and absolute phase 
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Returning to the eigen-values , we can write 

and t 9 with new values of a and (9 0 -9, ) as t,=me 
and t 2 = m tan 2 y > 1 1 


°1 

2jCPtP ) 
...( 2 . 2 . 16 ) 


whore is called the absolute phase, which dis- 
appears with power measurements. 


Thus the six parameters m,\^, ^ > f 3 completely 

determines the TSM.Eet us know express TSM m terms of 
these parameters or equivalently m terms of eigen-values and 
eigen-vectors . 


We can construct the unitary transformation - 


|"xl + xl x2 + xl 


— — 

1 0 

|xl + x2 x2 + x2_ 

— 

i. 


where orthonormal properties of xl and _x2 are used. 

Thus IF = (IF)"* 1 (oy) U" 1 = IF' = U + ,..(2.2.17) 


we can use the U to bring S into diagonal form 


U' SU = where Sfj =( ^1 21 


0 


...(2.2.18) 


Applying inverse transformation, noting that U 1 = U + * 


we get 


S = IF S d U T 


...(2.2.19) 


expressing U and S d m terms of eigen-values and eigen- 
vectors, we finally get. 


S= 


cosy- 

m 


sm* 


T 


m 


sm^ 


m 


cosy 


m 


' cos tn ~ sinT m 




*" jsint m CosT i 


m 


V 


me 


2jO+f) 


1 


0 


mtan ye 
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00S K ~ 

N sin K cos \ 


cos^_ 


m 


sm^ 

m 


L 


-smf cos**! 

m m 


...( 2 . 2 . 20 ) 


2.3 


Stokes Reflection Matrix 


Let a e.m. wave with stokes vector 

Pi ^I 1 , Q 1 , U 1 , v£| is incident on a target. 


Let the scattered wave have a stokes vector 
Ps £l s , Q S , U s y V^j . Then these two vectors are 
related by - 
■ si 


Ps = 


Q' 

'U fc 

V 


5 


mil 

mi 2 

ml 3 

mi 4 

m 2 l 

M2 2 

m2 3 

m24 

m3l 

m32 

m33 

m34 

[m4l 

m42 

m43 

m44 


i x 

Q 1 


u 


R Pi (2.3.1) 
= R Pi 


Where R is called stokes reflection matrix (also called 
Mueller matrix). 


Using law of reciprocity/ perrin C 1-4-]) showedthat the 

elements m = m . and hence the matax R is symmetrical. 

1 3 3 i 

R can further be rewritten in the following form- Huyren # 1970) 


Ao+Bo F C 

F -Ao+Bo G 

C G Ao+B 

H D -E 


H 

D 

-E 

Ao-B 


t" 


.. (2.3.2) 
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Rewriting that for a completely polarized wave, 
.2 2 2 2 
I = Q + U 1 +V 1 ' 

.2 2 2 2 

and for a partially polarized wave, I s > Q s +U S +v s * 
Huynen [ 1. ] arrives at the following conditions on the 

elements of R — matrix. 

2Ao (Bo+B ) ^ C 2 + D 2 


2Ao (Bo-B ) >, G z + 


2A0E 

2AoP 






CH - GD 


CG + DH 


(2 3 3) 


For a target, which scatters a completely polarized wave, 
the above relations satisfy with equality sign. In this, 
case, the target is called a single target and the corres- 
ponding stokes reflection matrix is denoted by M. On the 
other hand, if a target scatters a partially polarized wave, 
the above relations donot satisfy with equality sign. In 
this case the target is called a distributed target and 
corresponding stokes reflection matrix is denoted by R. 


For any arbitary targets, the stokes reflection 
matrix R associated ^ ith it can be decomposed into a single 
target M corresponding to completely polarized part and a 
noise target corresponding to unpolarized part. This theorem 
called target decomposition theorem, follows from the fact 
that a partially polarized wave with stokes vector (I,Q,U,V) 
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can be decomposed into a completely polarized wave with 

T 

stokes vector[l Q r U / v]and a unpolarized wave with stokes 

vector fl , 0,0,0] T where I = I + I and I i 2 = Q 2 + U 2 + V 2 
1 - u J c u c 

uniquely* 


The stokes reflection matrix R can be decomposed 
into M and M matrices as - 


— 

Ao+Bo 

F 

C 

H 


Ao+Bo 

T 

F 

C 

H 


Bo 

f n 

0 

0 

F 

-Ao+Bo G 

D 


T 

F 

-Ao+Bo G 

D 


f n 

B§ 

0 

0 

C 

G 

Ao+B 

-E ; 

““ 

C 

G 

Ao+B 

T -e t 

. + 

■ 0 

0 

b n 

-e n 

: H 

D 

-E 

Ao-B 


H 

D 

T 

-E 

Ao-B T 

1 


' 0 

0 

-e n 

-B* 






1 

— 



,1 






.. (2.3.4) 


where 


2Ao(Bo +B T ) 

CM 

Q 

+ 

CM 

O 

II 

Bo 


+Bo 

m m 

2Ao(Bo -B ) 

= g 2 + h 2 

B 

= b t 

T 

+b n 

N 

2AoE T 

= CH-DG 

. E 
and 

F 

= E 

T 

= F 

+E 

+f n 

2AoF T 

= CG+DH 





.^( 2 . 3 . 5 ) 


i ' 

The stokes reflection matrix M of a single target 
is irreduciabic m .he sense that M cannot bo further 
decomposed as sum of stokes reflection matrices corres- 
ponding to a collection ofsmgle targets. This therom 
points to a fundamental limitation of traditional attempts 
at " sectionalizmg" an object of complex shape into inde- 
pendent simpler shapes. Although at higher frequencies 
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when size to wave length ratio is large, such methods of 
computing R0S have had some success, the therom shows the 
futility of such attempts at lower frequencies. 

The target scattering matrix, lntEduced in the 
previous section and stokes reflection matrix describe the 
same target scattering characterstics . Hence one should 
able to get one matrix given the other one. The conversion 
follows from following steps. 

Voltage received at Rx terminals V = b'S a 
where a, b are Tx and Rx antenna polar izations, 
which converts voltage into electric field and 
vice-versa. 

* 

The received power is given by - P = V\f = 

s 

(B* Sa) (b* Sa)* ... <2.3.6) 

Simplification of above equation finally yields [ ] 

Ps = h (b) 1 M g (a ) 

where g (a), h(b) stokes vectors of Tx and Rx antenna 
polarizations , 

and M is the stokes reflection matrix defined above ..(2.3.7) 

The elements of stokes reflection matrix M is 
related to TSM elements as follows 



Let 


S 


— 

”1 
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x 1+j x 2 

Z 1 + J Z 2 

and M= 

mil 

ml 2 

ml 3 

ml 4 

x 1+: z 2 

Yi+jy 2 


m21 

m22 

m23 

m24 



I 

m3l 

m32 

m33 

m34 




m41 

m42 

m43 



mU = | [X 2 1+x 2 2+ y 2 1+ y 2 2+2z 2 1+ 2z^ ] 


m2 2 = 


m33 = ~ 


m44 = 


ml 3 = m3l =[^ 


Z 1 + Z 2 ~ X l y l “ x 2 y 2 

[x 2 + x 2 +y 2 +y 2 -2 z2.2 2 2 3 

2 2 

Z l +Z 2 + X l y l + X 2 y 2 

2 2 2 2 _ 

X 1 + X 2 - y l- y 2 ] 


...(2.3.9) 


ml 2 = m21 = 
ml 4 = m4l = 
m2 3 = m32 = 
m2 4 = m42 = 
m34 = m43 = 


z ± (y 2 - x 2 ) - z 2^ y l ” x l^ 

Z 1 ( x l~ ^2^ + Z 2 ^ X 2 + y 2^ 

z 2 (x x + y 1 ) p - z l ( x 2 +y 2 ) 
(x 2 y l - x l y 2 ! 

Z 2^ y 2 “ X 2^ + z l < ' y l' x i ^ 


Given the $tokes reflection matrix M, one can 
obtain the TSM uniquely exceot for absolute phase using 
the above relations. 
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Later when we compute scattering matrix for random 
rough surfaces, the elements of TSM are no longer deterministic 
but are statistically described. In these cases, the compu- 
tation of stokes reflection matrix can be done with slight 
modifications of above relations. Assuming that the 
scattering matrix elements are gaussian and neglecting the 
correlation bet’ scattering matrix elements, and bet 1 real 
and imaginary parts of each element, we can arrive at 
following relations for the elements of stokes reflection 
matrix R. 

mil = -| [(x 1 ) 2 +(x 2 ) 2 +(y 1 ) 2 +(y 2 )^+2(z 1 ) 2 +2(z 2 ) 2 +Var(x 1 ) 

+Var(x 2 ) +Var(y 1 ) + V-ar (y 2 ) \ 2Var <z_ L ) \2V&r( * 2 >3 

m2 2 = (z 1 ) 2 +(z 2 ) 2 - x 1 y 1 - x 2 y 2 + Var(z 1 ) + Var(z 2 ) 

m33 = ^['Z 1 ) 2 +(Z 2 ) 2 -Ay 1 ) 2 + ( y 2 ) 2 - 2(5 1 ) 2 - 2 d 2 ) 2 +Var (x 1 ) 

+Var (x 2 )+Var (y 1 )+Var (y 2 )-Var (z 1 )-2V ar (z 2 )] 
m 44 = ( z ) 2 t ( z 2 ) 2 t x^y 1 +x 2 y 2 +Var{z 1 ) t-Var (z 2 ) 

ml 3=m3l= \ [(^ ) 2 + ^ ) 2 - ( y 1 _) 2 - (y 2 ) 2 + VarCx^,) 

...(2.3.10) 

+ Var(x 2 ) - Var(y 1 ) - V a r(y 2 ) ] 
ml2=m21= z 1 (y 2 -x 2 ) - 5 2 C ^1"'*1 ) 
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ml 4 = 

m2 3 = 

m2 4 = 

m34 = 

where 


m4l = z 1 (x 1 +y 1 ) + i^+y,,) 
m3 2 = i 2 (it 1 +y 1 ) - ^(x^+y^ 
m42 = x 2 ^ - y 2 x x 

m43 = z 2 (y 2 -x 2 ) + 2i(y 1 ~x 1 ) 

bar denotes mean values . 
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CHAPTER 3 

SCATTERING FROM RANDOM ROUGH SURFACES 


In this chapter, we shall discuss the methods of 
computing the TSM and stokes reflection matrices associated 
with a target m the form of random rough surface, using 
the scattering theories. The expressions derived are 
extended to time varying surfaces. 


3 .1 Review of Scattering theories : 

T 

When an e.m. wave with stokes vector [ I 1 / q 1 , U 1 , V 1 ] 
is incident on a object, it induces electric and magnetic 
currents thro'out the whole volume of the body. These 
currents in turn produces a field away from the ^ ^ 

object. This field is called scattered fiEld 
with associated stokes vector [l s ,Q S , U S ,V^]^. 

In general, instead of taking into account 
currents in the whole volume of the object, 

one cap consider only the surface currents 

\ 

and compute the scattered field E ,H at f?^ 



some point. This assumption is perfectly valid m cases 
where conductivity is very large and penetration depth is 
small. The theories calculating scattered field using 
only surface fields are called surface scattering theories 


I 
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and these which calculate scattered field by considering 
fields in the whdb volume of the obj act is called volume 
scattering theories. 

The integral of Huygens-Kirchof f was modified 
to vector fields by Stratton-chu. The scattered field 
E at a point P away from the obstacle is expressed in 
terms of surface fields E and H as follows - 

^ S ( p) nxiti + (iv? J xvtp-Kft.E r ) vyy .(3.1.1] 

where ^ and h is the unit vector normal 

to the surface and E / H are fields existing on the surface 
of the scattering object and can be computed using appro- 
priate boundary conditions. The above equation is extended 
to Fraunhofer region by Slider (19 48). The expression for 
E S turned out to be the following — 

E S = BK 2 x vJ rSx£-V* 2,x(nx4)j exo^j v.?} c,‘ ...(3.1.2] 

where K 2 = ' K i = •*■/; ^ and 

^l / ^2 al " e wa-ve numbers m medium (l) & ( 2 ) 

and B = 1 k exp ( ikR) 

4 R 

The above equation/ when solved exactly gives 

.... - 

torn, h* * ^ 83^9151 
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accurate results. But nevertheless, one has to make 
some approximations m order to evaluate the integral. 

One such approximation, known as Kirchoff's 
approximation, essentially consists of approximating the 
boundary conditions. The approximation is as follows 

"If the radius of curvature at a point is large compared 
to wave length of the incident radiation, then the fields 
at that point can be approximated as the fields that would 
exist, if that point is a part of infinite tangent plane.* 

Using this approximation one can compute E and H 

— V cj 

fields existing over the surface and hence E can be 
computed. We shall use the above approximation to compute 

— -y S 

scattered field E . 


As the angle of incident increases, one has to 
take into account the shadowing effects. For our purposes, 
we shall use the shadow-function proposed by Smitht ] 
for Gaussian surf res. 


There is another method called perturbation method 
to compute the scattered field I s . This method is appli- 
cable only to surfaces with small-scale roughness. This 
method gives more accurate results than previous method 
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for slightly rough surfaces. This method is widely used 
to model the sea at HF/MF/LF. 

For some kind of surfaces, like sea at microwave 
frequencies, one has to use combination of above two methods. 
In this kind of surfaces, the small scale roughness is 
imposed upon large scale roughness. Hence one can compute 
the total scattered fields by computing the scattered field 
from small-scale rough surfaces and then averaging it 
over large-scale roughness. 


In our derivations, we will start with the 
scattered field expression derived by Bahar (1981) 


-a 

E S = Go j J Cor^T 1 U G?s ) exp ^ j On 3 -^ 1 ) * "rs? — d^'n ..(3.1.3) 


, _s 

r sh 


~LlH 

where E = 

E 

n 

E 



J E = 





j lV _ 


F 




, _HH, „W , _ , 

where F F are modified 

reflection coefficients. 


_HH €-1 

given by.- (cos 0 + cosy J2 

cos^G + s itL ) (l-^fe) 
i — 2 

(cos to + COS to-SH ) 


...(3.1.4) 
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where 0 is the angle of incidence 


is th 5 complex permittivity of the surface 
giver by 


€ r +j6o cr 


cos % - 1 , n = 


■' ''r _ lrv 


sm 

U(rs) = shadow function = 1 for illuminated area 

= 0 for non-illuminated 
area . 

1 £ 

T ,T are rotation matrices which rotates the 
incident field of reference-plane into local plane and 
vice-versa. They are given by - 



r— 




— 



T 1 = 

COS 

A 

- sm 

T f = 

cos 

0 

- sm 0 s 

i 

sin 

A 

cos 0 

1 ] 

7 

sxn 

0 
* s 

cos 0 s 


.. (3.1 .5) 


where 0 x is the angle bet’ local plane and incident plane 
and is the angle bet' 1 ocal plane and scatter plane. 

Assuming a mean s irface, m which E 1 is constant, one can 
make E 1 independent of area illuminated. Thus we can take 
E 1 out of the integral m equation (3.1.3*) and reqnte it 
as follows 

£ s = Go 0 jCoT^T 1 U(:r s ) exp jj(n S -n :L ) * r^ dA«n E 1 ..(3.1.6) 




wieow (rj(Q 
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comparing with E s = SE 1 , we can write, 

S = Goj j Q^T f FT 1 U(^ s ) exp £ j (t? s -£i ) . f s J dif.n ..(3.1.7) 

We will use equation (3.1.5) as the starting equation m 
ou' derivations in next section. 

3.2 Derviation of mean and variance of TSM 

for random rough surfacesT 

Let us consider a rough surface with roughness 
height varying m Z-directon and is denoted by Z=h(x,y)-» 

If h(x, y) is deterministic, the solutions are straight- 
forward and simple. But deterministic height variations 
are rarely realistic, model for practical applications. 

It is true that an individual area of terrain or seafurface 
at a give i instant of time, is given by a deterministic, but 
not a random function. But because of difficulty m deter- 
mining this function exactly and need to repeat the cal- 
culation for man/ other surfaces of same type, it is more 
convenient to rega d the considered surface as a particular 
realization of a random function with given statistical 
properties . 

In following derivations, we therefore restrict 
to the discussion entirely to random surfaces given by 
height distribution h(x,y). We make following assumption 


about the surface - 
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(1) The surface is assumed to be homogenous and 
isotnpic. Hence the correlation function 
C ^ x ( ^/Y c j) = < h(x,y) h (x+x^, y+y )> depends only 
upon displacement bet' the points where 
correlation is taken. As the surface is 
isotropic, the correlation function will only 

2 

depend upon the distance 1 = ( x^ + y^,) where 
X d'^d are displacement bet' the points m x 
and y - directions, regardless of its direction. 

(2) The co-ordinate system is defined in such a 
way as to make the x-y plane the mean of the 
surface, so that <h(x, y)> = 0. 

Let us start with the equation (3.1.3) given by - 

S= Go/JCoT^FT 1 exp |jk(i? S -rf 1 ' ) . "r^ U(r^ ) dJL*n ...(3.2.1) 

Let us evaluate mean and variance of S with following 
assumptions 

(1) The value of D = COT^FT 1 is constant through 

out the surface area illuminated. This assum- 
ption includes that the surface electrical 
properties are same thro' out the entire 
surface and angle# is taken as average angle 
0, which gives the average tilt of the 
entire surface. 
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(2) The function u(r s ), though dependent upon 

and y co-ordinates, ist taken as average 
value. 


Due to above two assumptions, D and U(r s ) 
can be taken out of the integral and we are left with 
exponential term to be integrated - 


S = Go D U(r ^) f j exp ^jk (n 3 -^ 1 ) . r^s^ dA-n 

t ~ 4 . v -ai \ v a* +v a* +v a* 

Let kln-n; = = xx yy z z 

= k(smQ cos^ s - smO cosA. ) 

s 

= k (smQ smi - sm8 smA ) 
s s 

= k (cos © s + cosG ) 

where (0 .A ) is the scattered direction and 


s ? s 


( © 


is the incident direction. 


..(3.2.2) 


For back scatter, 

© s = © , A s = > & = 

Hence, v x - 2k sm 0 

v 0 


y 

v_ = 2kcos© 

z 

where © is the incident angle. 

= x a x + y a y 


o 


+ h(x, y) 


.. .(3.2.3) 


a 


z 


The position vector r s 
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Thus k $!? - n 1 ) . r s = xv. + yv + h(x, y) v r 

x y ^ 


Thus , 


Co D U(r s ) 


X Y 

/ I exp (jv x +jv y+jh(x,y) z)dxdy 

x y 

-X -Y 


taking expected value, we get 


Go D X Y 


< s > = 


f f exp (jv x+jv y) < U(r ) exp ( jh (x, y )v, )> 


-X -Y 


x y 


For gaussian surfaces, 

< U(r_)exp jh(x,y)v > =KuC? )><exp(jh(x,y)v )> 
s z s z 


and identifying 


< exp(jh(x,y)v ) > = QC (v 


characterstic function, we get. 


< s > = Go D<U (r ) > 
mmm 1 s 


X Y 

(v ) / / exp(jv x x+jv y ) <3xdy 

Z — YV *• X 


, 0 v _ °o D<U (r )> 

^ ^ t smc (v X) smc (v Y) 'X (v ) 

A x y z 

Now let us calculate Var(s) s 


Var(s) , =<SS*> -<S><S^ 


< SS*> = Go D<(U(r )> 

-2 SI IS expAv x (x 1 -x 2 ) + jv y ( yi -y 2 )} 

^ -A, -j ' s. 

x <f exp jh 1 (x, y)-jh 2 (x,y) > dx 1 dy 1 dx 2 dy 2 

< S><S*> =|Go D< U (r )>| 2 , , \ / 

//exp jv x (x r x 2 ) + jy ( Yl -y 2 )C 

A 2 A 1 A ? _ k J 

dx 1 dy 1 dx 2 dy 2< 
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var(s) =< SS*> -<S><S*> 
= | Go D<U (r )> | 2 


ft $ ex K J VhW V«) ^ 2 


« -I 

~9i (v z ) X (v z ) dx 1 dy 1 dx 2 dy 2 
where x d = x 1 _x 2 ; y d = y 1 -y 2 . (v -v g )=<exp jh x (x, y)-h 2 (x t y) > 

let x d = lcosp and y d = Isxn p 


Substituting m above equation 
some simplifications 


, we obtain after 


var 


) = 


2 Go D < U(r )> | 2 co _ „ \ 

- 77 S - L ~t VW*** 1 ’.' * 


o xy * 

4C(v z ) x* (V z ) ] ldi 


For a gaussian surface, the height distribution is given by- 

W- (z) = exp (~z 2 /2 q 2 ) 

1 cfl[2% 

, N 1 , (a 1 -2c (l)z-, z 0 +z 0 ) 

W 2 z = 2nc^tl— c (l ) h “-i-J- 7-7“ 

2 cr ( 1 — c ( 1 ) ) 

where 'i ) is uhe correlation function. 

Therefore 


(v ) = exp (-%cr 2 Vo^ 

z z 

Ot 2 (v Z/ - V z ) = EXP ^ ® V 1 (1 "° (1>) ] 

Substituting m equation (3.2.4), we get 
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Substituting for v v andv from equation 

x / y z A 

(3.2.3) we get, 

< s>= Go D<U(r s ) ^ sinc (-2k sinSX) exp ( -2 d “k^cos^q) 

. .13.2.8) 

substituting (v 2 ), (v 2# -v z j values m equation (3.2.6) 

we obtain, 

2t j Go D<U(r )> | 2 « 

var(S) = * f J (v 1) f exp -or (1-c (l) ) 

A z J o xy L z 

A o 

-exp (-c 2 v 2 )] ldl *•(3.2.9) 

z 


Let us assume a gaussian correlation function. 

-1 2 /t 2 

Thus c(l) = e / 


where T is the correlation distance, 


2 2 

Let g = v z cr 


Thus exp f o" v_ (l-c(l))]= exp [-g(l-c Kl ) )] 

L Z **9 9 

co ^-mlVT 

expanding exp [-gc(l)] = £ 

m=o m 

—l 2 /t 2 

where we have substituted c(l) = e 


Thus the expression for variance becomes, 

o 2 /-p2 

2h|G D <b(r )> | ~ % _ a 00 

var(s) = — “ - J o (v xy l) e A 

A o 1 m=l 


mi 


ldl 


..(3.2.10) 


interchanging summation and integrations, we get. 
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7t|G 0 D <U(r s )> 


A' 


varfc) =— ^--2-^ e “ g Z £r exp ^-v^,TV4m 


m=l 


lH 

mm l 


2 _2 / * ) 
xy T /4m i 


... ( 3 . 2 . 11 ) 


The above equation is valid f or g «. 1 and g ^ 1, 

For the surface© with g>>l, instead of expanding 
exp {- gc(l)| , we expand the correlation function 

Thus c(l) = e- l2 / T2 CX. <1+1 2 /T 2 ) 


substituting in equation (3.2. 9 ) and simplification yields. 
2ft | G D <U(r )>| 2 


var(s) = 


A 2 g 


exp 


2 _2 
v xv T 


4v z & 


2 


(3.2.12) 


3.3 Shadow function. 


As we saw earlier , shadowing effects become very 
important, as we approach the grazing angles. In order to 
compute this function, we use the expression derived by 
Smith P for gaussian surfaces. In expressions 

derived m previous sections, we have taken shadowing into 
account thro' the function <U(r* )> 

Basically, <J(ig)> is the p^obablity that anoint 

not being shadowed. The shadowing m«3 udes the inviriblity 

of the point to the transmitter, receiver or both. The 

expression given by Smith is as follows — 

1 - 1/2 erfc ( 

“""ATpT'+T — 


u(F,e) = 
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2 

where \i = cot© W' = |C” (o)| and 

2 2 

2A(n) = (-|)l/2 Jet e-M- /2w _ erfc ( ll / V - 2w ) 

and 0 is the angle of incidence and c ^ (0) is the 
second derivative of the correlation function evaluated 
at origin. 


The above expression is valid for any correlation function 
, Let us evaluate the above function for a 
gaussian correlation function. 

0(1) . e- l2 A 2 

| c "( o ) | — 2/T 2 
and w = V" 2/T 


Substituting in above expression, we get, 

U(F,Q) = 1 ~ ** er£c 2 

A (cot© ) + 1 


«... ( 3 . 3 . 2 ) 


we will use this expressions in our equations derived 
above for < u (r^)> . 

< U(r s ) > = U(F, © ). 

Thus to evaluate the mean and variance of TSM, we 
should know following things^ 

1) physical properties of the surface (permittivity, 
permeability, conductivity); 
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2) anglo of incidence 

3) correlation function c(l) and correlation 
distance T 

4) surface roughness height deviation 

5) ^rea illuminated. 

3.4 Effect of temporal varial'on of the rough surface 
In the previous sections, we discussed scattering 
from a time invariant surface. In this section, we extend 
the results to include the time variation of the surface. 

In this case, the t_me variation o^ the surface is slow 
compared to the velocity of electromagnetic radiation. 

Hence we can apply quasi-static approximation to this case. 
By this approximation, all the equations derived in sec. 3. 2 
are still applicable with an additional term for doppler 
spectrtum. 

Let the height Deviation in z-direction is given by 
z = h(x,y,t«). In this case the equation for S(t) is 
modified as follows ° 

GoD xy 

S (t) = U(r s ^ J J exp (jv x x+jv y) exp jh( x,y , t) v z -dxdy 

X Y ..(3.4.1) 

and expected value of S(t) becomes, 

G oD 

<S(t)> = — t— %\v y j t) smc( v x)smc( v y) ...(3.4.2) 

** ^ x y 

whore ^c(v ,t) = bexp jh(x,y,t) v z > is the characterstic 

function, which is dependent on time. 
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The variance expression (3.2.6) gets modified time- 


varying surface as follows 

2 co 


va^S(t)J - jGoa^)>r_ 7 Jo(v i)Ettv z . - v 2 ,c(l),r)) 


XY 

-X(v z ,t] y fr (v z ,t)] ldl ...(3. 


4.3) 


where 1 = x“+y* and v^ y = v*+v^ 


The height deviation h(x,y,t) can be decomposed 
by using fourier integral as follows: 

h(x,y,t) =_j£/ dv x dv y h(v x ,v y ) e J t xv x +yv y- a ^ v x ,v y^ t ^ 

...03.4.4) 

The above representation indicates that the surface per- 
turations can be represented by monochromatic plane waves 
travelling with various wave vectors and frequencies 

Hence we obtain following representation for space- 
time correlation function: 


since c ( x^.y^T) = <h( x+xa) y +y d , t+t)h(x,y,t)> 
we get, 

oo 

e(^|.Y , d .i : < =/Jdv x c Jv y w(v x ,v y j exp[ j v x Hj+jv y y d 

-H(v,vJ ] ...(3.4.5) 

a, y 

where W(v ,v ) is a real, non-negetive function and is 
x y 

called rough surface spectrum. 

The frequency -ft- of each elementary harmonic is 
independent of the direction of its propagation for a 
homogenous and isotropic surface. It is related to the 
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spatial period 2 te/v , = (v ). The explicit form 

A / xy 

of this dependence (the ' dispersion- relation* ) is defined 
by the specific form of the equation of motion. 


Thus for a surface with dispersion relation (v ) 


the characterslic function becomes, 


, .. -1/2 cA, 2 sl/2 (vjt 

x(v_,t) = e ' 


v z e 


xy" 


2 2 , 


and X 2 (v z ,-v 2 ,c(1, )) = e~ a v 2 ( 1 " c ( 1 » )) 
Substituting in above equation (4.5 2), we get, 


var S(t) 


iGoD <U(r )>| 00 9 9 

b / Jo(v l)[exp -o v ( -c ( l , )) 

0 y 


A 


- exp( cr 2 v 2 ^ ldl 


Let us evaluate above equation for two cases. 


case (i) v a <<1, v a 1 


var S(t) 


| GoD<U(r )> 

— — o 


co 2 2 

f Jo(v xy l) e ff v z 


CL 


, 2 2\ m 0 
< v z g } c m (l, )-l] ldl 


m=o m 

GoD<U0r )>l^ 00 rf 2 2 2 2 m 

= ' 8 ' /Jo(v 1) e° v z ~ ( v z g ) c m (l, )ldl. 

o xy 


A 


Y 


m=l m 


interchanging sunniation 0 and integration signs, we get, 
,Gog<U(r s )>| 2 ,2,2 - (vVr 

var s(t) = ** — - 


-a v 
e z 


A 


m=l m 
.m/ 


Jj°(v l)c (]., )ldl. 

^ x y 


For m=l, this equation reduces to first-order perturvation 
solution Which is valid for v z <<1. For v z ° 1, we have to 


I 
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take sufficient number of terms and sum it to obtain 

var!js( t)^ . 

caso ( 11 ) v z <?>> 1. 

In this case instead of expanding the exponential, 
we expand the correlation function c(l,T). For the purposes 
of illustration, / r e shall consider first a 2-dimens lonal 
surface with z = h(x,t) and later extend the results to 
3~dimensional case. 


For two dimensional case, equation (3.4.3) is 
modified as follows i [ ] 

, . I GoD<U( r )>| 2 

varIS(t)U f e Jx d v x K 2 (v z ,-v z ,c(l,T.) 

i J A o 

- ^t(v z ,t) j£(v z ,t)] dx d ...(3.4.6) 


expanding the correlation function c(x^,X) near origin 
and taking only upto second order terms, we have, 

c(x d ,T) -- 1 + l / 2 c n *d + c 12*y +l/2c 22 2 

« 2 o(x ri -C) . d 2 c(x d , ) ...(3.4.7) 

where c n = ' c 12 “ " d x d dx 

d x d o 

^ 2 c(x d , ) 

C 22 = ' d '2 t ~ 

As v z a' i > 1, we shall neglect j^t(v z ,t)| term. 
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Thus equation (3.4. ) becomes, 

2 


var 




Gog< U(r]>| 




T 


/ e jv x x d v 02 2 


xp - v^ (l/2c,,x; 


il d 


+ c , 0 X 


12 x d +1/2 c 22 T") dx d ..(3.4.8) 


performing integration, by completing squares in the 
exponent and using the standard integral — 

2 


f e^ x dx = f % 


g- , vve finally obtain, 


var 


| s(t) j = 


| GoD<U( r ) > 


'T‘ 


A 


J2 2 
o v z 

exp 


fn 


2 ~ ~ 7 7 7UT 

v z {C ll J 


exp \ 




11 


>(c 


11 °22 


- c 2 ) 2 
C 12 ; 


...(3.4.9) 


since autocorrelation <s(t+ ) s*(t)> - var^s(t)^ for 
this case, the Fourier transform of above equation is given 
the spectral density. Thus, 

1 


s(to) a 




la - wj 

0xp ^ _ ~~-2“ 


(Aw)' 


2 2 


where (AwJ 2 = (c n c 22 - c 2 2 ) 


0 ) - v 12 

m x t — - 
C 11 


...(3.4.1o) 


The above equation is valid for m general, any 
dispersion relation S)~ . The scattered held spectrum is 


\ 
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a gaussian spectrum with mean and variance values dep- 
ending upon the surface roughness frequency, aspectangle 
and the kind of dispersion relation assumed. 


For three dimensional case, (i.e.) z 
the abO'/e equation modify to . - f ] 


S(to) = 


-n:|GoD<U(r )> 
kv z a |c n | 


exp 


1 

* Yntu 


h(x,y,t) 



where oa = ^1 

m “ 

and 

D= 



and (Au) 2 = ^ v z a ^ 

T - 


c xx 

c xy 

C 

xt 

c yx 

c yv 

c 

yc 

c X 

C 

y 

c z.x 

v x 

v y 

0 

c xx 

c xy 

C XT 

O 

X 

c yy 

C yt 


D 

d 


d = 


...( 3 . 4 . 11 ) 


XX 

**y 

yx 

c yy 


where the subscripts for C m indicates the corresponding 
second derivatives of c(xd,yd,X) where x d = y d =T = 


o 
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CHAPTER 4 

S“ “ 

APPLICATION TO PRACTICAL SURFACES 

In this chanter TSM and Stokes reflection matrices 
are numerically computed using the results obtained m 
sections 3.2 and 3.3, for the cases of an asphalt road 
and sea-surfaco. The Stokes reflection matrix obtained 
is decomposed using the decomposition thooren discussed 
m sec. 2. 3. The TSM corresponding to completely polarized 
component is calculated and the associated target parameters 
are obtained. The variation of these parameters .with 
frequency, aspect angle and surface roughness are studied. 
For sea-surface, the spectrum of the scattered field is 
also obtained. 

4. 1 Application to asphalt road ; 

The asphalt road has a complex permittivity of 
4.3 + jo. 15 and x-bancl and 2.5 + jo.65 at K£band[23]. 

The height dr ziation k a at those frequencies is assumed 
to bo 0.5 and tne correlation distance kT is assumed to 
be 0.5. Using tnese numerical values TSM and Smokes 
reflection matrices are computed as indicated above. 

The eigenvectors and eigenvalues of the problem S_x = tx* 
are determined from the solution Of the associated eigen- 
value problem S + Sx= x . This method of computing the 
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eigenvalues and eigenvectors is givenin. appendix A. 

From these, the associated target parameters are then 
computed and their variation with surface roughness and 
aspcctangle are plotted m the graphs [4.1.—] and [4.1*—]. 
As the linear polarization gives the maximum return, the 
values of XI and are constant and hence they all not 
plotted. 


4 . 2_ Application to sea-surface 

Sea surface, due to external forces like wind, is a 
time varying surface. Due to this, there is dispersion 
m frequency of scattered hold. 


The surface parameters of the sea are obtained from 
various sources as indicated below: 

Tne permittivity of sea-water is obtained utmg 
Jebye's formula given by - 


= € + 
oo 


£ 5 

■" 2Tir 

1+to > xr 


“■<(* 3-0 

i+tAr 


...(4.2.1) 


and^ r = 6 X - 3 


where € io 

O 



IS 

IS 


the static dielectric constant. 

the dielectric constant at infinite frequency 

the relaxation time 
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w = 2Tcf and f is the frequency at which 
dielectric constant is to be obtained 
a is the conductivity 
6 " 

o is tho free space permittivity. 

At 30°c, r = 7.?xl0 2 scc; £ s= 68.0 ,* ^ = ^.9J 
a = 1.567 mhos/m [20j 

The roughness height is taken to be the significant 
wave height. Significant wave height can be related to 

wind velocity as follows [19] 

-3 

H s = 7.15x10 W s 2 * 5 ...(4.2.2.) 

whore go is tho wind velocity m m/sec. 

o 

Tho correlation length T is estimated from the equation 

T = (0.556) 2 wj ...(4.2.3) 

Using above estimated values, the mean and variance 
of TSM and hence Stokes reflection matrix can be computed. 
The associated target parameters are obtained and their 
variations with aspect angle, frequency and wind velocity 
are obtained. It is found that m this case also, linear 
polarizations ga ro rise to maximum returns. Hence and 
are constants. Also it is observed that y= 0° for the 
sea-surface. Honce only the variations of m and y are 
plotted in Figs. 4.2 and 4.2.2. 


I 



52 


The temporal variations of the sea-surface give 

rise to the cloppler spectrum of the scattered waves. As 

discussed earlier, the dispersion relation ( v x » v ) 

the surface, modifies the surface spectrumW (v ,v ) of the 

x y 

sea as follows: [9] 

vV(v x ,v y ,T) = w ( v x » v y ) e“ jJl ^ v x ,v y^ T 

For back scatter, v = 2ksm© and v = 0. 

x y 


...(4.2.4) 


vtfe shall assume a gaussian spectrum for sea~surface. 


W ( v x »°) = p oxp (~v^T 2 ) ...(4.2.5) 

where T is the correlation length 

At low wind velocities when the wind is uniformly 
blowing the sea-surface can be assumed to have a linear 
dispersion. This assumption is valid, also for high wind 
velocities, at low frequencies. However, when the frequency 
becomes higher, or when the wind velocity is random, the 
dispersion becomes non—lincar. In order to study the 
spectrum for these two situations, let us consider- two 
kinds of dispersion relations, sin-e the surface is 
assumed to be homogenous and isotropic, dispersion relation 
depends only on v xy = liv x +v y and is equal to v x for back 
scatter case. 


case (l) linear dispersion • 

In this case, the dispersion relation is given by 


JL(v x ) = uv x = - 2kusin© where u 


is the velocity of the 
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ers 

scatter/ ,k is the wave number and 0 is the angle 
of incidence. 

ers 

The velocity of scatter/ is estimated from the 
relation given by ; [24] 

uU = fgv v + 0.04w c 1,5 + O.caV ...(4.2.6) 

X o s 

whore W the wind velocity is m/sec. and g is the acce- 
lOLation due to gravity. In the above equation, the first 
term corresponds to the effects cue to gravity, second 
term represents effect duo to orbital velocity and last 
term denotes wind effect. 

The spectrum is computed using the equations (3.1.10). 

It is found that, for this case, the spectrum is a gust a 

line shifted by the frequency o> m = t*v x « The variation of 

to with aspect angle, incident frequency and velocity of 
m 

the wind is plotted m graphs [4.2. — J to [4.2. — ] 

case (n) Non-lin ear di spersion % 

For non-linear dispersion, 

-A-(v x ) =uv x + av x 

where u = 

ers 

The value of y is taken as velodity of scatter/. 


and a = 


V x =° 


d 2 _l 

tv x=° 
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Since U-i has the dimensions of radians/sec. and has 
radians/m, a has a dimension of m /sec. In order to fit 
this dimension, the following relation for a is proposed: 

a = Dul' (4.2.7) 

where D is a constant p is the velocity of scatters and 
T is the correlation length. 

The velocity p is computed using the relation(4.2.6) . 
Using the above dispersion relation, the mean and variance 
of the scattered field frequency spectrum are obtained 
as follows ; 

w m = * v x 

and u m ) a =-§-(a^ ) (v z<J ) 2 (4.2.8) 

The standard deviation Am of the doppler is related 
to the half-power band width as [25]- Am m = 0.42 A co. .(4.2.9) 
where Am is the half-power bandwidth m radians/sec. 

The variations of the mean doppler w m and half-power 
bandwidth* w -/ith aspect angle, with wind velocity as 
parameter are plotted m fig. [4.2,3] to [4.2.8] 
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APPENDIX - A 

The eigen value problem Se = te* ....(A. I) 
can be solved as follows* 

Consider the conventional eigen vectors and eigen- 
values of the problem S + S_x = )\x ....(A. 2) 

since S is a complex, symmetric matrix, S + S is a Hormitian 
matrix. Hence, (1) The eigen-vectors of S + S matrix are 
orthogonal. (2) The eigen-values of S + S are real. 

Let us assume that ' S_x = e* ....(A. 3) 

Substituting (A. 3).n (A.2<), 

S + C* = S*0* =X X 

since S*o* = t*e, 

t*o = A x .... (A. 4) 

multiplying both sides by S, 
t*So = Sx 

1 1 1 = e* i 

A= 1 1 1 2 . ... (A. 5) - 

\ 

substituting m (A.4), we obtain, ,j 

t* e = 1 1 1 2 x 

i * Bu< [ 

e = t_x .....(A. 6) | 

Thus the relations (A. 3) and (A. 6) express, the 

j 

eigen-values of Che problem Se=te* om terms of conven- I 

+ 

tional eigen-values and eigen-vectors of the problem s 



CHAPTER 5 
CONCLUSION 


5,1 CONCLUSION; 

In this thesis an attempt has been made to model 
Radar- cl utter targets in terms of their polarization 
properties of the scattered returns. The scattered waves 
are decomposed into two parts - a completely coherent 
part corresponding to determinis tic part of the target and 
a incoherant part corresponding to a noise target. This 
deco reposition is carried out by Ihe use of Stokes Reflection 
Matrix associated with the target. The Target Scattering 
Matrix corresponding to the deterministic part of the 
target is obtained and a set of five phenomenological 
parameters associated with the target are also calculated. 

In order to illustrate above ideas, for a rough 
surface created as a random rough surface the mean and 
variance of Target Scattering Matrix associated it are 
computed. These values axe utilized to obtain Stokes 
Reflection Matrix. The expressions for mean and variance 
of the TSM are modified to include the temporal variations 
of the surface. The expressions derived are applied to an 
asphalt road and a sea-surface and the ’Five phenomenological 
parameters associated with these clutter- targets are obtamec. 
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The variation of these parameters with radar and surface 
parameters are also investigated. 

5.2 SUGGESTIONS FOR FURTHER WORK: 

In the above model, we evaluated the integral given 
m equation (3.1.7) making a number of approximations. A 
major approximation has been co lepiace the angles between 
local planes and incident pianos at different scattering 
points m the illuminated area by their average value. A 
more accurate and realistic approach would be to incorporate 
the angle-height dependence m cvaulatmg the integral. 

In co up u ting the Stokes Reflection Matrix, we have 
taken into account only the surface scattering component of 
the scattered field. This obviously restricts the applica- 
bility of the model to those sd tua Lions m which surface 
scattering is predominant or volume scattering effects are 
negligible. But clutter- targets like vegetation, sea-mce and 
show, where volume scattering effects are not negligible, 

A 

the Stolcos Reflection Matrix musl be computed by taking m 
to account the volume scattering effects [16,27]. 
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